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Voltage-gated proton (Hv) channels play an essential role in phagocytic cells by generating a hyperpolarizing 
proton current that electrically compensates for the depolarizing current generated by the NADPH oxidase dur- 
ing the respiratory burst, thereby ensuring a sustained production of reactive oxygen species by the NADPH oxi- 
dase in phagocytes to neutralize engulfed bacteria. Despite the importance of the voltage-dependent Hv current, 
it is at present unclear which residues in Hv channels are responsible for the voltage activation. Here we show 
that individual neutralizations of three charged residues in the fourth transmembrane domain, S4, all reduce 
the voltage dependence of activation. In addition, we show that the middle S4 charged residue moves from a 
position accessible from the cytosolic solution to a position accessible from the extracellular solution, suggesting 
that this residue moves across most of the membrane electric field during voltage activation of Hv channels. Our 
results show for the first time that the charge movement of these three S4 charges accounts for almost all of the 
measured gating charge in Hv channels. 



INTRODUCTION 

Voltage-gated proton (Hv) channels are found in phago- 
cytes, neurons, airway epithelia, muscle, and sperm 
(Decoursey, 2003; Okochi et al., 2009; Ramsey et al., 2009; 
Iovannisci et al., 2010; Lishko et al., 2010). In phagocytes, 
Hv channels play an essential role in bacterial clear- 
ing. Depolarization-activated H + currents through Hv 
channels in the phagosomal membrane compensate for 
NADPH oxidase currents, activated during a process 
called respiratory burst, ensuring sustained production 
of reactive oxygen species in phagocytes to neutralize 
engulfed bacteria. The importance of Hv channels is 
supported by the observation that inhibiting Hv chan- 
nels prevents further production of reactive oxygen 
species in the phagosomes. In addition, Hvl knockout 
mice are unable to efficiently clear bacterial infections, 
further showing the essential role Hv channels play in 
bacterial clearing (Ramsey et al., 2009). Because the volt- 
age dependence of Hv channels determines their physio- 
logical role, it is important to determine the molecular 
structure that defines the voltage sensitivity of Hv chan- 
nels. At present, it is unknown which charged residues 
underlie the voltage dependence of Hv channels. Here, 
using the limiting slope method and accessibility exper- 
iments, we have identified the charged residues respon- 
sible for voltage activation of Hv channels. 
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Each subunit in an Hv channel has four transmem- 
brane segments, called S1-S4 (Ramsey et al., 2006; Sasaki 
et al., 2006). These four transmembrane segments of 
the Hv channels are homologous to the first four trans- 
membrane segments in voltage-gated potassium (Kv) 
channels (Ramsey et al., 2006; Sasaki et al., 2006) . In Kv 
channels, the first four transmembrane segments form 
the voltage-sensing domain, whereas the fifth and the 
sixth transmembrane segments from all four Kv subunits 
together form the potassium-conducting pore domain 
(Jiang et al., 2003; Long et al., 2005). In contrast, Hv 
channels do not have a fifth or sixth transmembrane 
segment. In addition, unlike the tetrameric structure of 
Kv channels, Hv channels are composed of two subunits 
(Koch et al., 2008; Lee et al., 2008; Tombola et al., 2008). 
Deletion of the cytosolic domains removes the dimeriza- 
tion of the Hv subunits (Koch et al., 2008; Tombola et al., 
2008) . Interestingly, the monomelic Hv channel is a func- 
tional voltage-gated proton channel. Thus, each subunit 
comprises a voltage-gated H + permeation pathway (Koch 
et al., 2008; Tombola et al., 2008) . 

In Kv channels, the fourth transmembrane segment 
(S4) contains many positively charged residues and has 
been shown to function as the main voltage sensor in 
Kv channels (Aggarwal and MacKinnon, 1996; Larsson 
et al., 1996; Mannuzzu etal., 1996; Seoh et al., 1996; Yang 
et al., 1996; Tombola et al., 2006). S4 in Hv channels 
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contains just three charged arginine residues compared 
with, for example, seven charged residues in Shaker K + 
channels (Fig. 1 A) . A cysteine accessibility study of cyste- 
ines introduced in S4 of Hv channels shows that S4 un- 
dergoes a voltage-dependent transmembrane movement 
during gating, suggesting that the positive charges on S4 
are part of the voltage sensor that controls Hv channel 
gating (Gonzalez et al., 2010). However, that study only 
shows that S4 movement is correlated with channel open- 
ing and not that the S4 charges are the underlying gating 
charges. To demonstrate the role of the S4 positively 
charged residues with some direct methodology is im- 
portant given that other studies showed that individual 
neutralizations of the second and third S4 charges did 
not significantly decrease the voltage dependence of 
Hv channels, when measured as changes in the slope 
of conductance versus voltage curves at the midpoint 
of activation (Ramsey et al., 2006; Sasaki et al., 2006). 
Neutralizations of the first S4 charge increased the volt- 
age dependence of Hv channels in one study (Sasaki 
et al., 2006) but decreased the voltage dependence in the 
other study (Ramsey et al., 2006). In addition, Hv chan- 
nels truncated between the second and third charge in 
S4 still generated voltage-gated proton currents (Sakata 
et al., 2010). Therefore, the question of how much the 
three S4 arginines contribute to voltage sensing and ef- 
fective gating charge in Hv channels is open at present. 
Here we test the effect of neutralizing each of the three 
arginines in S4 to determine their contribution to the 
effective gating charge in Hv channels by the more accu- 
rate method of measuring the gating charge in the limit 
of very small channel open probability (i.e., the limiting 
slope method; Aimers, 1978; Sigg and Bezanilla, 1997). 

MATERIALS AND METHODS 

Mutagenesis and expression 

Site-directed mutagenesis, cRNA transcription, and injection of 
cRNA encoding Ciona intestinalis Ci-VSOP (here called Ci-H v ) 
into Xenopus laevis oocytes were performed as described pre- 
viously (Sasaki et al., 2006; Bruening-Wright et al., 2007). The 
ANAC Ci-H v has a stop codon at Val270 and initiator methionine 
replacing Glul29 (Gonzalez et al., 2010). 

Cysteine accessibility 

Cysteine accessibility to MTS reagents in excised inside-out 
patch clamp and two-electrode voltage clamp (TEVC) experiments 
was assessed as described earlier (Gonzalez et al., 2010) using 
MTSEA (2-aminoethyl MTS), MTSET (2-(trimethylammonium)ethyl 
MTS), or MTSPT (3-(trimethylammonium)propyl MTS). TEVC 
external solution contained (in mM) 88 NaCl, 1 KC1, 1 MgCl 2 , 
1 CaCl 2 , and 100 HEPES, pH 7.4. Currents were recorded using a 
Dagan GA-1B amplifier (Dagan Corporation), low-pass filtered at 
1 kHz, and sampled at 5 kHz. Oocytes were injected with 50 nl of 
1M HEPES, pH 7.0, to minimize pH changes caused by the pro- 
ton currents. This results in a concentration of ^T00 mM HEPES 
in the cytosol. 1M-HEPES injection prevents the internal pH 
from changing significantly by the proton currents (Gonzalez 
et al., 2010). Both pipette and bath solutions for excised patch 



recordings contained (in mM) 100 HEPES, 2 MgCl 2 , and 1 EGTA, 
pH 7 with NMDG Currents were recorded using an Axopatch 
200B amplifier (Molecular Devices), low-pass filtered at 1 kHz, 
and sampled at 5 kHz. To increase the proton currents in some 
recordings, we used pH, 5.5 (and Mes instead of HEPES). The 
rate of MTS modification was only determined using patches 
where solution exchange was <100 ms (e.g., Fig. 4 D; Larsson 
et al., 1996). The rate of modification was measured by plotting 
the change in the current by the MTS reagent as a function of the 
exposure to the MTS reagents (exposure = concentration x time, 
measured in [mM s]; e.g., Fig. 4 G). The current versus exposure 
plot was then fit with an exponential to get the second order rate 
constant (unit/mM/s) of the MTS reaction. 

Limiting slope theory 

For a two-state voltage-gated channel (one closed state and one 
open state), the gating charge z can be estimated by fitting the 
conductance versus voltage, G(V), curve by a Boltzmann: G(V) = 
G 0 /(l + exp( — z(V — Vi/ 2 )/kT)). Alternatively, one can plot the 
G(V) curve in a log-lin plot and fit the slope at small open prob- 
abilities, i.e., the limiting slope technique (Aimers, 1978; Sigg and 
Bezanilla, 1997). This is because at low open probabilities, log 
G(V) — > log G 0 + zV/kT. For channels with more states than two, 
the fit of the G(V) with a Boltzmann only gives a lower estimate of 
the total gating charge, z total , needed to move to activate the chan- 
nel. For example, a Hodgkin-Huxley type channel with only two 
subunits and a gating charge z = 2 eo per subunit (Fig. 1 B, inset) 
has a G(V) = G 0 /(l + exp(z(V - V 1/2 )/kT)) 2 . The G(V) from this 
two-subunit channel is not well fit with a Boltzmann (Fig. 1 B, red 
dashed line) , and the fit gives a lower z value (z flt = 2.43 eo) than 
the real total gating charge per channel Zt oB i (z^ = 2 e 0 + 2 e 0 = 4 e 0 ) . 
Fitting the G(V) with a Boltzmann raised to a power of two (z fit = 
2 e 0 ) or three (z flt = 1.87 e 0 ) fits the G(V) much better (Fig. 1 B, 
solid red and dashed teal line, respectively). Fitting the G(V) with 
a Boltzmann raised to a power of two or three is equivalent to as- 
suming that the channel has two or three identical subunits (in a 
Hodgkin-Huxley type channel), respectively. In this case, the fitted 
z fit is per subunit and the total estimated charge per channel z total 
is two times z flt for the Boltzmann to the second power (z total = 2x2 
eo = 4 e 0 ) and three times z flt for the Boltzmann to the third power 
( z totai = 3 x 1.87 e 0 = 5.6 e 0 ). From this it is clear that the fits of the 
G(V) with a simple Boltzmann underestimate the true total gating 
charge, whereas fitting with higher powers of a Boltzmann could 
lead to overestimates of the true gating charge. In contrast, the 
limiting slope technique is less model dependent and will for most 
models give a good estimate of the equivalent gating charge, if 
the slope is measured at low enough open probabilities. For ex- 
ample, for the two-subunit model, the log G(V) — > log G 0 + 2zV/kT 
for low open probabilities. Using the limiting slope technique for 
simulated data down to a open probability of 10~ 3 , the fit to the 
three-state channel gives a z Umit = 3.8 e 0 (Fig. 1 C, teal line), close 
to the correct value, z total = 4 e 0 . 

Limiting slope experiments 

To measure the limiting slope, currents were recorded from 
excised inside-out patches containing WT or charge-neutralized 
Ci-Hvl channels in response to slow ramps of voltage (1-2 mV/s) 
from a holding potential of —60 to 40 mV (or 0 mV) and back to 
— 60 mV (e.g., Fig. 1 D). Solutions for excised patch recordings 
contained (in mM) 100 HEPES, 2 MgCl 2) and 1 EGTA, pH 7 with 
NMDG, in pipette and 100 Mes, 2 MgCl 2 , and 1 mM EGTA, 
pH 5.5 with NMDG, in bath solution. Currents were recorded using 
an Axopatch 200B amplifier, low-pass filtered at 1 kHz, and sam- 
pled at 2 kHz. Up to 20 current traces were averaged and low-pass 
filtered at 5 Hz to improve the signal to noise. There were no 
significant endogenous currents in response to the voltage ramps 
(Fig. 1 D, teal line), most likely because of the almost complete 
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Figure 1 . Limiting slope estimates the effective gating charge. 
(A) Alignment of the S4 region (teal box) of Hv channels from 
C. intestinalis (Ci-Hv), Xenopus (XenHv), Mus musculus mouse 
(mHv), Homo sapiens (hHv), Callus gallus (ChHv), and Danio rerio 
(DanHv) and S4 from Kvl.2, Shaker K, and domain IV of Navl. 4 
channels. The arrow shows residues R1454 in Navl. 4 and R368 
in Shaker that have been shown to move completely across the 
membrane during channel activation (Yang et al., 1996; Baker 
et al., 1998). Red letters indicate positively charged residues in 
S4. (B) Simulated conductance versus voltage, G(V), curve (black 
circles) for a three-state, two-subunit Hodgkin-Huxley type model 
with a gating charge of z = 2 e 0 per subunit. The G(V) was fit with 
a Boltzmann curve, G(V) = G0/(1 + exp(-z(V - V 1/2 )/kT)) (z flt = 
2.43 eo; dashed red line) , and a Boltzmann curve to the second 
(zfi, = 2.0 e 0 , resulting in a z tota] (fit) = 2 x z fit = 4.0 e 0 ; red line) and 
third power (z fit = 1.87 e 0 , resulting in a z toIal (fit) = 3 x z flt = 5.6 e 0 ; 
teal dashed line). (C) Simulated data and fits from B shown in 
lin-log plot (colored as in B). Also the best fit to the limiting slope 
at negative voltages is shown (log G(V) = constant + z limit V/kT; 
ziimit(nt) = 3.8 e 0 ; teal straight line). (D) Currents from excised 
patches containing WT Ci-Hvl channels (black line) or from 
patches taken from uninjected oocytes (teal line) in response 
to slow voltage ramps (1.5 mV/s) from —60 to 0 mV and then 
back to —60 mV. (E) The calculated conductance (I/(V — E reT )) 
during the up and down ramp (data from D). Also shown is the 
calculated conductance from a patch from an uninjected oocyte 



absence of any ions permeable to endogenous oocyte channels in 
these patch solutions. Currents were leak subtracted for ramps 
between —100 and —80 mV. We calculated conductance versus 
voltage, G(V), curves by dividing the currents by (V — E rev ). The 
conductance-voltage curves from up and down ramps were indis- 
tinguishable, as expected when the channels are near gating equi- 
librium during the measurements (Fig. 1 E). G(V)s of WT Hvl 
channels were not well fit by a simple Boltzmann (Fig. 1, F and G, 
dashed red lines). Instead, we plotted the conductance on a semi- 
logarithmic scale (Fig. 1 G), and the G(V)s were fitted in different 
voltage intervals to identify the limiting logarithmic slope, as de- 
scribed earlier (Gonzalez et al., 2010). 



RESULTS 

The effective gating charge is reduced 
by S4 neutralizations 

We studied the Hvl channel from the sea squirt C. intesti- 
nalis, Ci-VSOP or Ci-Hvl , because these channels can be 
expressed at higher densities in Xenopus oocytes than 
their mammalian Hvl homologues. The higher densities 
of Ci-Hvl channels expressed in oocytes enabled us to 
more reliably measure the proton currents at low open 
probabilities, as necessary for the limiting slope tech- 
nique. The effective gating charge of WT Ci-Hvl chan- 
nels is 5.9 eo (dimer) and 2.7 eo (monomer) , as measured 
from the limiting slope conductance (see Materials and 
methods; Gonzalez et al., 2010). The effective gating 
charge is the amount of charge that needs to move rela- 
tive to the transmembrane electric field to turn on the 
proton conductance. The higher value in the dimeric 
channels is consistent with cooperative activation of the 
proton conductance in the two subunits (Gonzalez etal., 
2010; Tombola et al., 2010). Here, we replaced, one by 
one, the three arginines on the S4 segment of Hv chan- 
nels from C. intestinalis with asparagines (Fig. 1 A) to de- 
termine the contribution of each of the three arginines 
present in S4 to the gating charge in Hv channels. 

To measure the effective gating charge in the argi- 
nine mutants, we applied very slow voltage ramps (1-2 
mV/s) to inside-out patches from oocytes expressing 



(open circles) . The currents were leak subtracted in response to 
ramps between —100 and —80 mV. The reversal potential E rev was 
estimated to —85 mV in this experiment (pH = 5.5 and pH D = 7). 
Note that the conductance estimated during the up ramp over- 
laps with the conductance during the down ramp, showing that 
the channels are at equilibrium at all voltages. This shows that the 
speed of the ramp is slow enough to correctly estimate the conduc- 
tance. (F and G) Conductance, G, measured as in D from excised 
patches containing WT Ci-Hvl channels in response to slow volt- 
age ramps (1.5 mV/s) from —60 to 50 mV. The G(V) was fit with 
a Boltzmann curve, G(V) = G0/(1 + exp(-z(V - V V2 )/kT)) (z flt = 
3.5 eo; dashed red line), and a Boltzmann curve to the second 
power, G(V) = G0/(1 + exp(-z(V - V 1/2 )/kT)) 2 (z fit = 3.0 e„, re- 
sulting in a z total (fit) = 2 x z fit = 6.0 e 0 ; solid red line) . (G) Data and 
fits from F shown in lin-log plot (colored as in F). Also the best 
fit to the limiting slope at negative voltages is shown (log G(V) = 
constant + z limit V/kT; teal straight line). z limit (fit) = 5.3 e 0 . 
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the different S4 charge neutralization mutants and 
measured the resulting currents (Fig. 2 A) . The current- 
voltage curves from up and down ramps were indistin- 
guishable (Fig. 2 B), as expected when the channels are 
near gating equilibrium during the measurements. We 
calculated the conductance from these current mea- 
surements and plotted the conductance versus voltage 
curves in semilogarithmic plots (Fig. 2, C-E). The limit- 
ing slope at negative voltages has been shown to be an 
accurate measure of the effective gating charge neces- 
sary to open the channel (Aimers, 1978). A substantial 
decrease of the effective gating charge as measured 
from the limiting slope conductance was measured when 
R255, R258, or R261 was neutralized compared with the 
effective gating charge in WT Ci-Hv channels (Fig. 2 F) . 
Double and triple neutralizations of the three S4 charges 
did not produce functional channels or their expres- 
sion was too low to conduct patch clamp experiments. 
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Figure 2. Neutralization of any S4 charge reduces the effec- 
tive gating charge. (A) Currents from excised patches containing 
R255N Ci-Hvl channels in response to slow voltage ramps (2 mV/s) 
from —60 to 40 mV and then back to —60 mV. (B) The current 
during the up and down ramp (data from A). (C-E) Conduc- 
tance (G) calculated as G = 1/(V — E rev ) from measurements 
as in A and B from excised patches in response to slow voltage 
ramps (2 mV/s) for R255N (C), R258N (D), and R261N (E) Hv 
channels in pH; = 5.5 and pH„ = 7. The limiting slope fitted as 
G = G 0 / e zSVF/RT is shown for each recording. (F) Effective gat- 
ing charge (z8) estimated from the limiting slope for charge- 
neutralized and WT Hv channels (n = 4) from experiments as in 
A and B. Error bars are SEM. 



The reductions in effective gating charge for the sin- 
gle neutralizations suggest that each one of three S4 
charges contributes significantly to the voltage sensing 
of Hv channels. However, a decrease in the effective 
gating charge caused by the charge neutralizations 
could also be the result of a change in the dimerization 
probability or extent of cooperativity of the two Hv sub- 
units. Therefore, we also made the S4 charge neutral- 
ization mutations in the form of Hv channels that had 
the cytosolic domains deleted (ANAC; Koch et al., 
2008), a deletion which prevents dimerization (and 
thereby also abolishes cooperativity) . In these mono- 
melic Hv channels, the effective gating charge decreases 
substantially when R255, R258, or R261 is neutralized 
(Fig. 3), as if all three S4 charges contribute to the ef- 
fective gating charge in monomeric as well as dimeric 
Hv channels. 

The S4 movement is curtailed by S4 neutralizations 

The reduction in gating charge produced by the neutral- 
izations is surprisingly large (>2 eo for dimeric channels 
and >1 eo for monomeric channels), especially for the 
R255N mutation (Fig. 3 C) . A reduction of 1 eo per sub- 
unit is the maximally expected reduction in the gating 
charge by the neutralization of one charged residue, as- 
suming that the charge neutralization only removes a 
charge from S4 and does not alter anything else in the 
channel. If one adds up the gating charge reduction for 
all mutations, one gets a total of 12 eo, a number which is 
twice the 6 eo measured for WT Hv channels. This result 
suggests that the effect of the neutralizations is not a sim- 
ple removal of the charged residue. If, in addition, the 
charge neutralization also reduces the distance moved by 
the remaining charges in the transmembrane electrical 
field, then the contribution of the remaining charges to 
the gating charge will also be reduced. We therefore 
tested whether some charge reduction in Hv channels 
also could be caused by an altered S4 movement, espe- 
cially for the R255N mutation. 

We used cysteine accessibility to test for altered S4 
movement in the R255N mutation. There are two na- 
tive cysteines in the N terminus and one native cysteine 
in the C terminus of C. intestinalis Hvl. However, we 
have previously shown that the currents from WT Ci-Hvl 
are not modified by the application of internal or ex- 
ternal MTS thiol reagents, such as MTSET or MTSES 
(2-sulfonatoethyl MTS; Gonzalez et al., 2010). There- 
fore, we used the WT C. intestinalis Hvl channel as the 
background for our cysteine mutations and cysteine ac- 
cessibility experiments. I262C has earlier been shown to 
be accessible to MTSET from the intracellular solution 
in closed channels (Fig. 4, A and B; Gonzalez et al., 
2010). Modification of I262C by intracellular MTSET 
abolished the current through I262C channels (Fig. 4 B) . 
We combined I262C with R255N and tested the inter- 
nal accessibility of I262C in the presence of R255N 
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(Fig. 4, C-G) . Using a fast perfusion system (Fig. 4 D) , 
we were able to apply 1 mM MTSET for 1 s to the in- 
tracellular side of the patch. This allowed us to apply 
MTSET specifically while the channels were held at 
a specific voltage. We monitored the current from 
R255N/I262C channels in response to a 100-mV pulse 
every 10 s, while applying 1 mM MTSET to the cytosolic 
face of an excised patch for 1 s every 10 s at —100 mV 
(Fig. 4 E) . In contrast to I262C channels (Fig. 4, A and B) , 
internal MTSET did not modify closed R255N/I262C 
channels (Fig. 4, C-G) . Even applying internal MTSET 
at such a negative voltage as —140 mV did not modify 
R255N/I262C channels (Fig. 4 G, triangles) , as if S4 
does not move as far inward toward the cytosol in 
R255N/I262C channels compared with I262C channels 
(Fig. 4 H) . The decreased range of S4 motion is consis- 
tent with the large reduction of the effective gating 
charge caused by the R255N neutralization (Fig. 2 F). 
A possible reason for the decreased range of S4 motion 
is that the two end points of S4 movement are determined 
by charge-charge interactions of S4 charges with the 
two clusters of negative residues at the intracellular and 
extracellular border of the channel (Fig. 4 H; Long 
et al., 2007). When the most extracellular S4 charge is 
neutralized (R255N) , S4 would only travel as far inward 
as to make the next S4 charge (R258) interact with the 
internal negative cluster (Fig. 4 H) . I262C would in this 
S4 position be located at the same place as residue 259 
in WT channels, which we previously showed to be inac- 
cessible from the intracellular solution (Gonzalez et al., 
2010) . A similar mechanism with a reduced outward S4 
motion could explain the large reduction in gating 
charge for the R261N neutralization (Fig. 2 F). 

The middle S4 charge moves completely 
across the membrane 

A simple interpretation of the reduction of gating charge 
for the S4 charge neutralizations (Fig. 2) is that these S4 



charges move across most (if not all) of the electrical field 
during channel activation (Fig. 5, A and B), suggesting 
that these S4 charges move completely from a cytosolic 
accessible position to an extracellular accessible position. 
However, the accessibility experiment above shows that 
charge neutralizations might alter the extent of the S4 
movement, thereby complicating the interpretation of the 
charge neutralizations. To get further independent evi- 
dence that these S4 charges move from the cytosolic side 
to the extracellular side of the membrane, we also tested 
the accessibility of R258C, the middle S4 charged residue, 
to internal membrane-impermeable cysteine reagents 
(MTSET) in closed and open channels. R258 lies in be- 
tween two residues, L256 and V259, which we earlier have 
shown to be inaccessible from both the intracellular and 
extracellular solutions in both closed and open channels 
(Gonzalez etal, 2010). However, mapping these residues 
(L256, R258, and V259) onto existing models for Kv chan- 
nels in closed and open states suggests that R258 could be 
accessible from both sides of the membrane, even if 256 
and 259 are not accessible from either side of the mem- 
brane (Fig. 5, C-F) . We therefore decided to test the ac- 
cessibility of R258C from both sides of the membrane. 
External accessibility of 258C was tested by applying MTS 
reagents on whole oocytes expressing R258C Hvl chan- 
nels, whereas internal accessibility of 258C was tested by 
applying MTS reagents on excised inside-out patches from 
oocytes expressing R258C Hvl channels. 

Internal MTSET modifies R258C channels (Fig. 6, 
A and B). MTSET slowed die time course of the R258C 
currents (Fig. 6, A-C). Internal MTSET modified R258C 
channels when applied at negative potentials (k = 0.63 ± 
0.04/ mM/s; n = 7) but not at positive potentials (Fig. 6, 
D and E) . In contrast, external MTSET did not modify 
R258C. However, externally applied MTSEA, which is 
smaller than MTSET, modified R258C channels (Fig. 7, 
A and B). The membrane-impermeable MTSPT, which 
has a longer linker between the charge and the thiol 
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Figure 3. Neutralization of S4 charges reduces the effective gating charge in monomeric Hv channels. (A) The current during 
up and down voltage ramps for monomeric R255N ANAC Ci-Hvl channels. (B) The calculated conductance (I/(V — E rev )) during 
up and down voltage ramps for monomeric R258N ANAC Ci-Hvl channels. The conductance is similar during the up and down 
ramps, showing that the ramp speed is slow enough to estimate the steady-state conductance. (C) Effective gating charge (zS) es- 
timated from the limiting slope for charge-neutralized and WT Hv channels in dimeric (red) or monomeric (teal) form (ANAC; 
see Materials and methods; n = 4). Error bars are SEM. 



Gonzalez et al. 279 



A 



C 0.02 
g 

O 0.00 



c 



Before MTSET 





I 












pZZ — 











B 

I262C 



E 0.02 



u 0.00 



After MTSET 



0 500 1000 1500 2000 0 500 1000 1500 2000 



Time (ms) 



D 



Time (ms) 




R255N/I262C 

0.02, MTSET 



o 



500 1000 1500 2000 

Time (ms) 





1000 2000 3000 4000 

Time (ms) 



After MTSET 



0 1000 2000 3000 4000 

Time (ms) 

G H 
• « ■ i 



500 1000 1500 2000 

Time (ms) 
I 

R255N/I262C 




Exposure (mM s) 



262C- 

100 mV +60 mV -100 mV +60 mV 



Figure 4. R255N alters the extent of S4 movement. (A and B) 
Currents from an excised inside-out patch containing I262C 
channels in response to voltage steps from — 40 to 80 mV, before 
(A) and after (B) application of internal MTSET. 1 mM MTSET 
was applied for 10 s on closed channels held at —100 mV in be- 
tween the two recordings. (C and F) Currents from an excised 
inside-out patch containing R255N/I262C channels in response 
to voltage steps from —40 to 100 mV (in 20-mV increments), 
before (C) and after (F) 1 mM MTSET was applied 10 times 
at — 100 mV for 1 s. pH; = 5.5 and pH 0 = 7. (D) Currents in response 
to a 100-mV voltage step every 10 s during MTSET application on 
open R2355N/I262C channels. The speed of the perfusion of the 
cytosolic face of the patch was tested by applying 1 mM MTSET 
(in pH 7.0 solution) at 100 mV for 1 s, followed by washout at 
100 mV (pH 5.5 solution). Wash-in and washout were monitored 
by the fast changes in current amplitude induced by changes in 
internal pH: pH; = 5.5 in rinse and pH ; = 7.0 in MTSET solution. 
pH D = 7.0 (pipette solution). (E) Currents in response to a 100-mV 
voltage step every 10 s during MTSET application on closed 
R255N/I262C channels. Closed-state modification was tested by 
applying 1 mM MTSET at — 60 mV for 1 s, followed by washout. 
Same patch as D. (G) Normalized currents in response to voltage 
steps to 100 mV after sequential 1-s applications of 1 mM MTSET 



reactive end than MTSET, also modified R258C (Fig. 7, 
C and D). External MTSEA modified 258C channels 
when applied at depolarized potentials (k = 0.18 ± 
0.03/mM/s; n = 3) but not at hyperpolarized potentials 
(Fig. 7 E). The external accessibility data suggest that 
R258C is located in a narrow extracellular crevice in the 
open state. The accessibility of R258C from both sides 
of the membrane, but with opposite voltage depen- 
dence, shows that R258C moves from a position acces- 
sible from the cytosolic solution in closed channels to a 
position accessible from the extracellular solution in 
open channels. 

DISCUSSION 

We have shown here that neutralizations of individual S4 
charges reduce the effective gating charge in Hv chan- 
nels, as if all three S4 charges contribute to the gating 
charge. In addition, we show that during Hv channel 
opening, the displacement of the middle S4 charged resi- 
due (R258) encompasses the whole membrane electric 
field. When channels are closed, this residue is accessible 
from the cytosolic solution. In contrast, when channels 
are open, this residue is accessible from the extracellular 
solution. Assuming that the membrane electrical field 
falls over the inaccessible portion of the protein, this sug- 
gests that R258 contributes 1 eo per subunit to the gating 
charge. In voltage-activated Na channels, it has earlier 
been shown that a cysteine introduced at the homolo- 
gous residue to R258 (Fig. 1 A, arrow) can bind MTSET 
on the extracellular side of the membrane and then the 
intracellular reducing agent TCEP (Tris (2-carboxyethyl) 
phosphine) could remove MTSET on the intracellular 
side, clearly showing that S4 residue modified with MTS 
reagents can move from one side of the membrane to the 
other side of the membrane (Yang et al., 1996). Our 
present data, together with our earlier work in which we 
showed that the movement of S4 precedes Hv channel 
opening (Gonzalez et al., 2010), strongly support a 
model in which the voltage dependence of Hv channels 
is caused by the movement of the three S4 charges. 



applied at -100 mV on I262C (red boxes) and R255N/I262C 
(black boxes) channels, or sequential 1-s applications of 2 mM 
MTSET applied at -140 mV on R255N/I262C (black triangles) 
channels. MTSET was only applied during the inter-episode time 
at a holding voltage of —100 mV (or —140 mV; as in E). Error 
bars are SEM. (H and I) Cartoon showing proposed interactions 
of the positively charged S4 residues (R255, R258, and R261) with 
the two clusters of negatively charged residues in S1-S3 at the in- 
tracellular and extracellular borders of the Hv subunit in I262C (H) 
and R255N/I262C (I) channels. We propose that it is the inter- 
actions of the S4 charges with the two negatively charged clusters 
that set the limits for the extent of S4 movement (Fig. 3 [A and B] 
from Gonzalez et al., 2010). 
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Earlier studies of the roles of the three S4 charges in 
the voltage dependence of Hv channels were inconclu- 
sive or contradictory (Ramsey et al., 2006; Sasaki et al., 
2006) . The conclusions of these studies were based on 
changes in the Boltzmann fits of the whole conductance 
versus voltage G(V) curves from different charge-neutral- 
izing mutations (Ramsey et al., 2006; Sasaki et al., 2006). 
However, the number of gating charges estimated from 
Boltzmann fits of G(V) curves are only lower estimates of 
the true gating charge. This is because mutations could 
change the slope at the midpoint of the G(V) curves by 
altering the cooperativity between channel subunits and/ 
or altering the relationship of the different activation 
steps in a multistep activation pathway of the channels, 
and not just by decreasing the number of gating charges 
(Zagotta et al., 1994). Therefore, it is very unreliable to 



use changes in Boltzmann fits to estimate the contribution 
of individual charges. Instead, here we used the more reli- 
able method of measuring the limiting slope at low con- 
ductance to estimate the effective gating charge for the 
different mutations (Aimers, 1978; Sigg and Bezanilla, 
1997). However, even the limiting slope technique has 
limitations. One of the possible errors is to not have mea- 
sured the conductance at low enough open probability 
to have reached the true limiting slope. For example, in 
Navl.4 channels, the effective gating charge was esti- 
mated to 4—7 eo per channel from limiting slope experi- 
ments on macroscopic currents where the conductance 
could be estimated down to an open probability of 10 3 
(Oxford, 1981; Stimers etal, 1985), whereas the effective 
gating charge was estimated to 12 eo per channel from lim- 
iting slope experiments on single channel experiments 




Figure 5. Cysteine accessibility is consistent with an S4 movement in which all three S4 charges significantly contribute to the effective 
gating charge. (A and B) Proposed S4 movement and charge-charge interactions in WT Hv channels according to present and pre- 
vious (Gonzalez et al., 2010) cysteine accessibility data mapped onto the voltage-sensing domains from the Kvl. 2-2.1 chimera channel 
structure (Long et al., 2007; A) and a closed-state model of Shaker K channels (Pathak et al., 2007; B). Positively charged S4 residues 
(Rl, R2, and R3) and negatively charged glutamate residues (El and E2) are shown as ball and stick (except R2 = R258, which is space 
filled and colored yellow) . Purple residues are only accessible from the extracellular solution in open channels and blue residues from 
the intracellular solution in closed channels (Gonzalez et al., 2010). Red residues (256 and 259) are not internally or externally accessible in 
closed or open channels (Gonzalez et al., 2010). Solid lines indicate proposed lipid bilayer boundaries, and dashed lines indicate 
proposed MTS accessibility as the result of water-filled crevices. Proposed positions of the charges of the three S4 arginine residues rela- 
tive to the inaccessible portion of S4 are indicated. Assuming that the electrical field falls linearly over the inaccessible portion of S4, 
R255, R258, and R261 would contribute 0.80 e 0 , 1.0 e 0 , and 0.80 e 0 , respectively, to a total gating charge Q= 2.6 e 0 . This is similar to 
the contributions associated with these charges in FEP/MD simulations on the Kvl. 2 voltage-sensing domain (Khalili-Araghi et al., 
2010). (C and D) Two side views of the S1-S4 of an Hvl channel monomer in the open (C) and closed (D) state rotated 180 degrees 
relative each other are shown in space fill. (E) Top view of the open state shown in C. (F) View from intracellular solution of the closed 
state shown in D. R258 is shown in yellow, and residues 256 and 259 are shown in red. R258 is accessible from the extracellular surface 
through a water-filled crevice (C and E) and from the intracellular surface through a water-filled crevice (D and F). Neither residue 
256 nor 259 is visible from either side of the membrane in either closed or open states, consistent with earlier accessibility experiments 
(Gonzalez et al., 2010). Arrows show aqueous access to R258. 
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where the conductance could be estimated down to an 
open probability of 1CT 7 (Hirschberg et al., 1995). In 
principle, there is no way to know that the true limiting 
slope has been reached in any particular experiment, un- 
less there is another way to estimate the gating charge 
movement (e.g., from gating currents) . Our experiments 
are limited by the signal to noise at very low open 
probabilities for macroscopic currents because of the 
very small single channel conductance of Hvl channels 
(Cherny et al., 2003) and because no gating currents 
have been measured from Hvl channels. However, the 
estimated effective gating charge of 6 eo in the Ci-Hvl 
channel from our limiting slope experiments matches 
well with the estimated charge transfer caused by the S4 
movement measured by our cysteine accessibility experi- 
ments (Gonzalez etal., 2010). 

The neutralizations of the S4 charges caused surpris- 
ingly large reduction in the effective gating charge for the 
charge neutralization mutants (Fig. 2) . A reduction in the 
measured effective gating charge of only 1 eo per subunit 
is expected for any charge neutralization if the neutraliza- 
tion does not cause any other change to the protein or the 
applied transmembrane electrical field. One possibility is 



that we have not reached the true limiting slope for these 
mutations because of the limiting range of open probabili- 
ties in our experiments for these mutations (Figs. 2 and 3) . 
However, in that same open probability range, the slope 
of the G(V) curve for WT Ci-Hvl channels has already 
reached a value that is within 10% of the limiting slope 
value measured at even lower open probabilities (Fig. 1, 
E and G; Gonzalez et al., 2010), suggesting that we have 
not grossly underestimated the gating charge for the mu- 
tants. Another possibility is that the assumptions that the 
neutralization does not cause any other change to the pro- 
tein or the applied transmembrane electrical field are not 
necessarily true for all types of mutations. One should be 
especially concerned about the validity of these assump- 
tions for mutations of charges that are essential for the 
conformational changes underlying channel gating. If 
one adds up the gating charge reduction for all mutations, 
one gets a total of 12 eo, a number which is twice the 6 eo 
measured for WT Hv channels. However, similar neutral- 
ization experiments on Shaker K channels also reduce the 
total gating charge by >1 eo per subunit for some muta- 
tions (Aggarwal and MacKinnon, 1996; Seoh etal., 1996). 
If one adds up the gating charge reduction for all Shaker 




Figure 6. Internal MTSET modifies R258C channels only in the closed state. (A and B) Currents from an excised inside-out patch con- 
taining 258C channels in response to voltage steps from —60 to 60 mV, before (A) and after (B) application of internal MTSET (1 mM) 
for 10 s and washout for 10 s. The holding potential was —60 mV. (C) Currents in response to a voltage step to 40 mV before and after 
application of internal MTSET. The current time course was fitted with an exponential (red lines) with t = 70 ms before MTSET and 220 ms 
after MTSET. (D) Currents from an excised inside-out patch containing R258C channels in response to sequential voltage steps to 
40, 100, and -100 mV. V holding = -80 mV. pET = 5.5 and pH„ = 7. 1 mM MTSET was applied for 1 s during the -100-mV voltage step. This 
protocol was repeated every 10 s for 10 times (only four traces are shown for clarity) . The inset shows currents during the 40-mV voltage 
step, showing the slowing of the current time course in response to consecutive MTSET applications. (E) The rate of MTSET modifica- 
tion at —100 mV (closed boxes) or 100 mV (open circles) was measured by fitting the current time course during the 40-mV voltage 
step (from data like in panel D) to a double exponential with t1 = 70 ms (unmodified channels) and t2 = 220 ms (modified channels). 
The amplitude of the slower component is plotted versus the cumulative MTSET exposure. The amplitude of the slower component is 
assumed to be proportional to the number of MTSET-modified 258C channels. The data were fitted with an exponential to yield the 
second order rate constant k for the MTSET reaction with R258C. k = 630/M/ s for MTSET applied at — 100 mV. No significant modifica- 
tion was seen when MTSET was applied only at 100 mV. n = 7. Error bars are SEM. 
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mutations, one gets a total of 24 eo, about twice the 12-14 
e 0 measured for WT Shaker K channels. This suggests that 
the effect of the neutralizations is not a simple removal 
of the charged residue. If, in addition, the charge neu- 
tralization also reduces the distance moved by the re- 
maining charges in the transmembrane electrical field, 
then the contribution of the remaining charges to the 
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Figure 7. External MTS reagents modify R258C channels only 
in the open state. (A and B) Currents from an oocyte express- 
ing 258C channels in response to voltage steps from —50 to 50 mV, 
before (A) and after (B) application of external MTSEA (150 uM). 
The holding potential was — 60 mV. MTSEA was applied at 40 mV 
for 65 s in between the two recordings, showing that 258C is 
accessible to the external solution in open channels. pH D = 7.4. 
(C and D) Currents from an oocyte expressing 258C channels in 
response to voltage steps from —50 to 50 mV, before (C) and after 

(D) external application of membrane-impermeable MTSPT. The 
holding potential was —60 mV. 1 mM MTSPT was applied at 40 mV 
for 65 s in between the two recordings, showing that 258C is 
accessible to the external solution in open channels. pH c = 7.4. 

(E) Sequential 5-s applications of 150 uM MTSEA followed by 10-s 
washout was applied at the indicated voltage, and then the current 
was monitored in response to a 40-mV step every 30 s. The rate of 
MTSEA modification at —80 mV (closed boxes) or 40 mV (open 
boxes) was measured using the time constant for the current activa- 
tion (left y axis) for the first 1,000 ms after the start of the 40-mV 
voltage step. The activation time constant versus the cumulative 
MTSEA exposure was plotted and fitted to an exponential. As a 
control, 150 uM MTSEA was applied at 40 mV on WT Ci-Hvl chan- 
nels, and the time constant for activation was monitored (open tri- 
angles; right y axis) . No change in the activation time constant was 
noticed for WT Ci-Hvl channels. Error bars are SEM. 



gating charge will also be reduced. Indeed, in Shaker K 
channels it has been shown that some S4 charge neutral- 
izations alter the extent of S4 movement (Baker et al., 
1998) . Also, in recent extended molecular dynamics (MD) 
simulations of Kv channels, it was shown that an S4 charge 
neutralization reduces the travel of S4 (Jensen et al., 2012) . 
This could explain why the reduction of the total gat- 
ing charge could be >1 eo per subunit (Aggarwal and 
MacKinnon, 1996; Seoh et al., 1996) . We therefore tested 
whether some charge reduction in Hv channels also could 
be caused by an altered S4 movement, especially for the 
R255N mutation. Indeed, we show that the extent of mo- 
tion is smaller in R255N channels compared with R255 
channels. The curtailed movement of S4 in channels with 
the R255N mutation is most likely a contributing factor to 
the large reduction in effective gating charge for R255N 
because the smaller S4 movement leads to a smaller con- 
tribution of the two remaining S4 charges to the effective 
gating charge. 

In Fig. 5, we propose a model for S4 movement in 
Hv channels that is consistent with the accessibility data 
and the limiting slope measurements of the charge neu- 
tralizations. We base this model on the accessibility data 
and available crystal structures and molecular models 
of related Kv channels (Long et al., 2007; Pathak et al., 
2007; Gonzalez et al., 2010). The contribution of each 
S4 charge can be estimated by assuming that the voltage 
is falling linearly over the inaccessible part of S4 (Fig. 5) . 
This assumption results in an estimated total contribu- 
tion of 2.6 eo for the three S4 charges (Fig. 5). A more 
realistic estimation of the contribution of each S4 
charge for a similar model of Kvl.2 channels has been cal- 
culated using free energy perturbation (FEP)/MD 
simulations without the assumption of a linear voltage 
drop (Khalili-Araghi et al., 2010). In the Kvl.2 model, 
R2 contributes 0.7 eo, R3 contributes 1.0 eo, and R4 
contributes 0.8 eo to the total gating charge in an iso- 
lated voltage-sensing domain. In our model of Hv 
channels, R1-R3 (R255, R258, and R261) correspond 
to R2-R4 in the Kvl.2 model. Assuming that Hv and 
the isolated voltage-sensing domain of Kvl.2 share 
similar structures, this model generates a gating charge 
per subunit of 2.5 eo. Both the simple estimation based 
on a linear voltage drop and the more realistic FEP/MD 
simulation give a total gating charge (2.6 and 2.5 eo, re- 
spectively) from the three S4 charges that is in reason- 
able agreement with the measured gating charge per 
subunit of 2.7 eo in Hv channels (Gonzalez et al., 2010), 
suggesting that no other charges contribute signifi- 
cantly to the gating charge in Hv channels. In conclu- 
sion, our data show for the first time that the three S4 
charges are the main voltage-sensing residues respon- 
sible for the voltage dependence of Hv channels and 
that the extent of movement of these S4 charges is suf- 
ficient to explain the voltage dependence of opening 
in Hv channels. 
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